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Abstract 

The Hydrodynamics of Superfluid Turbulence (HST) describes the flows (or counter- 
flows) of Hell in the presence of a chaotic set of vortex filaments, so called superfluid 
turbulence. The HST equations govern both a slow variation of the hydrodynamic vari- 
ables due to dissipation related to the vortex tangle and fast processes of the first and 
second sound propagation. This circumstance prevents an effective numerical simulations 
of the problems of unsteady heat transfer in Hell. By virtue of a pertinent multi-scale 
perturbation analysis we show how one can eliminate the fast processes to derive the 
evolution equation for the slow processes only. We then demonstrate that the long-term 
evolution of a transient heat load of moderate intensity obeys the nonlinear heat conduc- 
tivity equation, often referred to as the Dresner equation. We also compare our approach 
against the Dresner phenomenological derivation and establish a range of validity of the 
latter. 
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Introduction and scientific background 

As it is known, upon exceeding some critical value of the counterflow velocity in case of the 
motion induced by the heat load applied to one end of a channel, a chaotic set (or tangle) of 
quantized vortex filaments emerges in superfluid helium. The presence of the vortex tangle dra- 
matically changes all the hydrodynamic properties of Hell in comparison with those prescribed 
by usual two-fluid hydrodynamics (See for details book by Donnelly [T] and review article by 
Tough 12]). Due to a small value of the critical velocity such situation takes place almost in any 
experiment on studying hydrodynamic and heat transfer phenomena in Hell. Consequently, 
any flow or counterflow of HeH should be investigated with the help the set of hydrodynamics 
equations coupled with Vinen equation for the vortex line density dynamics (see |2|) . This set 



1 



of equations known as the Hydrodynamics of Superfluid Turbulence (HST) have been obtained 
in papers of Nemirovskii and Lebedev Geurst [S] and Yamada et. al. within the scope 
of different approaches. The Hydrodynamics of Superfluid Turbulence as well as the methods 
for its investigation are outlined in the review paper by Nemirovskii and Fiszdon [Zj. The 
set of HST equations is extremelly cumbersome, therefore, there is no wonder that to achieve 
quantitative results in virtually any case of practical interest one is bound to turn to numerical 
methods. Yet, numerical simulation of nonstationary flows of Hell faces one serious obstacle. 
The point here is that the set of HST equations that unifies the hydrodynamics equations of 
superfluid helium and the Vinen equation for the vortex line density evolution is initially of 
hyperbolic type. As a result, a slow variation of hydrodynamic variables due to dissipation 
is accompanied by the fast processes related to propagation (and possible reflections from the 
boundaries of the channel) of the first and second sounds. If one is interested in the slow evolu- 
tion of the temperature, velocities and the vortex line density only, particular details of sound 
propagation and manifold reflections come out completely irrelevant, yet requiring rather ex- 
tensive numerical resources. Indeed, a typical time of transient processes is of order of seconds, 
whereas a typical mesh size in numerical investigations of acoustic phenomena should be of 
order of at least microseconds to account for the forming of the shock fronts. Thus, it seems 
attractive to try to get rid of the fast modes by a pure analytical procedure. In the present paper 
we realize effective separation of the slow from the fast modes using multi-scale (in the case 
multi-time) asymptotic perturbation techniques (See, for example, book by Nayfeh |H]. The 
main idea here is to introduce several time scales, each of which responsible for the description 
of a corresponding process. As an example we recall, perhaps, the simplest possible case of the 
overdamped harmonic oscillator. Given a certain relation between the frequency and the damp- 
ing decrement, sinusoidal oscillations die away shortly, whereas slow exponential relaxation may 
take quite a while to die out. Therefore, the long-term dynamic behavior changes drastically. 
The same is true for our case: dimensional analysis of the HST equations shows that there is a 
dimensionless criterion, Strouhal number Sh that is the ratio of the counterflow decrement to 
inverse time that takes the heat pulse to cross the channel. Depending on the input parameters 
of the problem the number Sh can be either large or small. The latter circumstance allows us 
to come up with an effective procedure for separation of the fast from the slow modes. The 
paper is organized as follows. The second section is devoted to dimensionless analysis of the set 
of HST equations. In the third section we carry out the multi-scale asymptotic perturbation 
procedure for the case of heat load of moderate intensity, of order of a few W/cm'^ i.e. for 
Sh >> 1. We then establish that the slow processes are governed by the nonlinear heat con- 
ductivity equation similar to that derived phenomenologically by Dresner P]. We also discuss 
the criteria of validity of the Dresner equation and the range of its applicability. In conclusion 
we resume the obtained results. 

Dimensionless equations of HST 

As it was noted above the full set of the HST equations is very complicated and in general, 
does not permit any analytical investigation at all. Having in mind to study processes related 
to heat transfer in Hell, we propose the following simplifications. 

(i) First, we restrict ourselves to the case of counterflow of Hell which occurs whenever one 
applies a heat load to a reservoir filled with superfluid helium. 

(ii) Second, we study only quasi-one- dimensional cases, i.e. either pure one- dimensional (d=l). 
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or cylindrical (d=2), or spherical (d=3) geometries. 

(iii) Next, we neglect by the nonlinear convective terms on the left-hand sides of the equations 
of motion. Usually these terms come in play in the cases of very large heat fluxes (when no 
separation of the "fast" and "slow" mode is possible at all) and are responsible for the forming 
of the shock fronts. 

(iv) Finally, we drop the dissipation function which is very small in the cases to be studied. The 
above assumptions correspond to the situations brought about in the overwhelming majority 
of experiments on transient heat transfer. To describe the counterfiow of Hell it is convenient 
to choose the following set of variables: dimensionless velocity of the normal component , 
dimensionless temperature T' and dimensionless square root of the vortex line density (VLD) 



G {G = \J C/ Coo 1 where Coo is equilibrium with respect to the normal velocity value of VLD). 
Under the assumptions listed above the set of equation of HST (see for details e.g. review 
article of Nemirovskii and Fiszdon [Zj) is reduced to the following form: 

with the following dimensionless variables: 

t = ^t', X = Lx', Vn = VnO V'n. T = ^T', C = CooG' . (4) 
U axU 

Here L is a length of the channel, U is the velocity of the second sound, Vno is a characteristic 
value of the normal velocity at the boundary Vno = q/ pcT (here g is a heat flux, which is 
actually a function of time, a is entropy per unit mass and ar is a derivative of entropy with 
respect to the temperature). 

At this stage it is worth making a few remarks concerning the set of equations (^-(jSl)- If 
one used dimensional variables Vn,T and C instead of their dimensionless counterparts, one 
would easily recognize the equations widely used in studying the dynamics of intense heat 
pulses. Nevertheless a somewhat awkward use of variable G in lieu of the vortex line density is 
necessary to develop the multi-scale perturbation procedure. 

It is easy to see from set ©-© that the dynamics of the hydrodynamic variables is specified 
by several dimensionless criteria, namely the Vinen parameter Vi, the Strouhal parameter Sh 
and the Mach number Ml. They are defined as: 

Sh = -— , Ml = -f , Vt = — - 5 

UTd U A{T)Pspn P 

Here we have introduced the following notations: a and f3, the parameters of the Vinen 
equation; A{T) the Gorter-Mellink constant; Vl is a drift velocity of the vortex tangle; r^, a 
characteristic time of attenuation of the second sound due to the vortex tangle introduced by 
Vinen (see jS]): 
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Numerical estimation of parameters (0) shows that within the temperature interval 1.4 - 
2.1 K parameter Vi changes from 0.2 to 1.5. It does not depend on applied heat flux q (or on 
quantity Vno )■ The Mach number, although dependent on the heat flux, remains very small. 
As far as parameter Sh is concerned, it strongly depends on the heat flux and can be either 
large or small compared to unity: 

Sh << 1, for small heat flux (7) 

Sh » 1, for large heat flux (8) 

Thus, the main parameter, drastically affecting the behavior of the system, is the Strouhal 
parameter Sh which is nothing else than the ratio of the counterflow damping decrement due 
to an interaction with the vortex tangle to the inverse time of the flight of the heat pulse. 
Conditions ((Zj)-( IH} can be used for effective separation of the processes with characteristic 
time to ~ L/U (second sound behavior) from those with characteristic time ti ~ {L/U) /Sh 
(evolution due to counterflow- vortex tangle mutual friction). 

Having in mind to investigate the cases of practical interest, we estimate Sh numbers for 
heat load of moderate intensity q = 1 -r- 10 W/cm"^ and channel sizes L = 10 -r- lO^cm. 
Using the definitions (jSj), (jH)) of the Strouhal number as well as thermodynamical parameters 
we obtain that in the temperature region T = 1.4 -h 2. IK the Strouhal number is of order of 
10^ — 10^, i.e. much greater than unity. Therefore, we can restrict ourselves to studying the 
cases of large Sh numbers. 

Multi-time scale method 

In this section we carry out elimination of the fast processes. As mentioned above we will 
concentrate on the case of large Strouhal parameter Sh » 1. Put in another way, attenuation 
of the second sound is assumed so strong that the convective regime of heat transfer quickly 
gives way to another regime, reminiscent of nonlinear heat conduction in usual newtonian fluids. 
Formally this leads to the degeneration of the initially hyperbolic equations into a parabolic 
equation describing the slow evolution of the velocity and temperature fields. To demonstrate 
this transition explicitly we invoke a pertinent multi-scale asymptotic perturbation theory (see 
e.g.Nayfeh [H]). Following this method, we introduce different time scales. 

to = t'; t[ = et'; t'^ = eH' , (9) 

where e = 1/Sh « 1. We look for a solution to the set of the HST equations (HJ-Q in the 
form of an asymptotic series 

K = \/o'(a;',^o,^'i,^2) + eK(^',^o,i'i,4) + e'^'(a:',to,t'i,t2) + 

T' = T^(x',t[„t;,t2/) + eTi'(a;',t;„t'i,t^) + e2T^(x',t^,t;,t^) + .... (10) 

G = Go{x', tp, t'^, t2) + eGi{x', tg, t'^, t2) + e^G2{x' , tp, t'^, tg) + .... 

In accordance with multi-scale perturbation analysis (see (H]) the coefficients in the series for 
dimensionless normal velocity V^, dimensionless temperature T' and qauntity G are supposed 
to have an order of smallness 0(e°), i.e. to be of order of unity. The simple chain rule follows 
from relations - ((EI): 

d d d ^ d . . 



The next step in study of the slow evolution of the heat pulse consists in substituting the 
multi-time scale series (Pl-lfTUj) into equations of HST (P)-©- Gathering terms of the same 
order of magnitude with respect to e we come up with a chain of equations leading to divergent 
(secular) solutions. Canceling step by step these secularities we then obtain a hierarchy of 
equations of different orders in parameter e, governing different stages of the evolution of the 
fields of temperature, counterflow velocity and the square root of the vortex line density. 

Let us develop the outlined scheme in detail. Excluding the temperature variable T' from 
the set of equations (IH)-(ini) we rewrite it in the as: 



V' - 



dt' y " 




2 (dG dG\ 



"^-1 dx' "/ 



(12) 



(13) 



Implementation of the above procedure depends on a particular statement of the problem 
under consideration. To be specific, we assume that at the moment a heater is switched on 
the thermal front starts out propagating into undisturbed bulk of helium which is, in fact, 
the case in most experiments. Such statement, however, brings up some apparent paradox. 
Indeed, substituting series ()1U|) for G and V and using rule (fTT|) we find out that to zero order 
in e, quantity G is equal to lo'(a;', tg, t'^, ^2)- Put in another way, it follows that VLD C{t) 
immediately takes its equilibrium value with respect to the normal velocity. But that is wrong 
as there actually is some finite time Ty of development of the VLD (see [S]). Hence, we seem to 
end up in quite a contradictory situation. To find a way out, let us compare the characteristic 
time of the vortex tangle development Ty against the characteristic time of propagation of heat 
pulse L/U. As it is well known (see j3] ), the former is determinated by relation: 



a{T) q 



-3/2 



Taking into account a scatter of data it follows 

g3/2^ ~ 40 ^ 140 



1^3/2 



(14) 



(15) 



It is seen from (15) that for heat pulses of moderate amplitude 1 \W/cm'^] < < 10 \W/cni?] 
propagating in channels of length L ~ 10 -r lO^cm , which are typical of experiments on 
transient heat transfer in Hell, the quantities Ty and L/U are of the same order of magnitude. 
It implies that for the time of order (L/U), typical of the slow processes, the vortex line density 
has enough time to adjust to the velocity field, i.e. to take its equilibrium (with respect to the 
relative velocity) value C^o = ^^(^)^- In other words, if we start our scheme not at t = 
when there are no vortices in the {)ulk of helium but at time t = Ty when the VLD has adjusted 
to its equilibrium value, we can then proceed without any discrepancies. At times smaller than 
Ty when G is small relations ()12j ) -f [T^ can be rewritten as 



'_d_ 
dx' X 



Id-l 



d_ 
dx' 



X 



id-l 



dt'^ 



v 







(16) 



5 



The meaning of relations (fTB|) and (fT7|) is obvious. For times smaller than L/U the heat 
pulse propagates according to wave-like equation because of the absence of vortices ( excluding 
extremely small background value). The vortex line density, being very small in comparison 
with its equilibrium value jCoo, changes mainly by first generating term of the Vinen equation^. 

Let us now proceed with the general scheme. The following step in studying the slow 
evolution of the heat pulse consists in substituting multi-scale variables ITUj) into equations 
( fT ^ - (fT^ . Gathering terms of the same order, we wind up with the chain of equations: 

To order e°. 

Go = K. ^^^^ 



To order e^. 



V - r — 'Ml. 1 ^ 

— Vi dx> 

^o^^(K + 2Gi) = -^Gl + [-^^{x"^'^ V^)) . 



Using (fTHj) the second equation of set (fT^ can be integrated over t'^ to yield: 

(20) 



•^^0^1 — J\-^ ^^1^^2) + Vi dx' + 



did (^td~l 9_^3 

dx' x"^-^ dx'y-^ ^0) at'i 



It is seen that for I'q — > 00 the quantity VI — > 00 as well. In order to cancel this divergence 
we have to require the quantity inside the brackets to be equal to zero. This procedure is called 
the cancellation of divergencies in secular terms in an asymptotic series. Having accomplished 
this procedure, we figure out that to the first order in the evolution of the heat pulse obeys 
the set of equations below: 

K = v/(x',t;,4), (21) 



2Ml l_ dV^ 

'^0 



Thus we have derived a set of equations which describes slow variations of hydrodynamic 
variables ( with the values of parameters discussed above). It can then be inferred from (jl8p 
and (|23|) that the leading part of the normal velocity depends only on the slow time and is 
governed by the parabolic-type evolution equation. This is the main result of the paper. It 
holds true provided that all the necessary creteria for application of the outlined perturbation 
procedure are met (see discussion above). Later we shall establish the range of validity of our 

^To see what would exactly happen at times smaller than ry one had to resolve a fairly difhcult problem of 
propagation of the heat pulse which creates the vortices on its way and interacts with the very vortices. The 
corresponding investigations are described in the review article |7| 
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approach. At this point, however, we complete the calculation. To this end let us study the 
next, order. We follow exactly the same procedure as above, cancelling divergences one step 
at a time. Leaving out simple but somewhat tedious calculations we write down the final result: 
To order e^ . 

VI = viixyp 

Gi = Gi(x',t]^) 

Go = Vq .24) 
V' - r - 1 avo' ^ ' 

•^1 ~ Vi dx' 

_d_\//3 _ did ( \ri\ 

dt', ^0 ~ dx' x"i-^ dx'y-^ 



dt'^ '0 dx' dx' 

q _9_n/'2T/') — 9 1 d (^ld-1 T//\ I AMl 1 SV^ 

dt'-^^^Q ^l) ~ dx' x"'-^ dx' ^1) ^ Vi dx' 

It is plain that the last equation accounts for a small correction to the velocity field and 
hence is of little interest to us. 

The Dresner nonlinear heat conductivity equation 

Let us consider the first order in greater detail. With the goal to compare our results 
to the well-known Dresner heat conductivity equation(see e.g. [H]) we consider here only the 
simplest geometry {d = 1). Combining equation [21] with the results of the previous section 
we arrive at the set of dimensionless equations for the normal velocity, and the temperature: 

Here V^',T', are zero terms in series (fTUI) for the dimensional normal velocity and the tempera- 
ture. We recall that to this order in e the vortex line density takes its equilibrium value with 

/^'^. Next, excluding quantity 



respect to normal velocity , e.i. C = V^^ . Next, excluding quantity from the set of equations 
written above we obtain the relation: 

dt[ dx' \dx' J ^ ^ 

Relation ()27|) formally coincides with the widely used nonlinear heat-conductivity equation 
derived by Dresner In this connection it is worth discussing the method used by Dresner. 
He started with the Gorter-Mellink relation which in our notation transforms to: 

3T' 

— = -Sh V'l (28) 

Recalling that heat fiux q is related to normal velocity Vn via q = STVn, he concluded that 
the heat fiux was proportional to the cube root of the temperature gradient. 



q oc 



dx' 
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Furthermore, using the energy conservation law, 



dt dx 

Dresner derived an equation similar to equation ()27p. His approach, however, is not self- 
consistent from point of view of the full set the HST equations. Indeed the Gorter-Mellink 
relation (|28l ) corresponds to a steady-state and its direct use for a nonstationary case is not 
valid. The correct equation should be obtained from equations (1), in which quantity G is 
replaced by its equilibrium value V'n- 

dt' ^ dx' ~ ^^^"^ 

Throwing out this term implies that = V^{x) which, in turn, means that the temperature 
is also a function of x only and does not change in time. Accordingly T' connected to via 
Gorter-Mellink relation does not change in time. 

In our approach dVQ/dt'^ ^ and there is no contradiction. Our careful consideration 
explicitly demonstrates that the regime of the nonlinear heat-conductivity equation takes place 
only for "slow" time (~ L/U) and validity of the overall procedure requires meeting the criteria 
dwelt upon above. The proposed allows also to determine a range of validity of the Dresner 
equation (j2Zj). One restriction is that slow regime takes place at a time scale of order {L/U) 
therefore a parabolic type heat-transfer equation works at times greater than t^. The latter 
then provides the lower limit of the appplicability of the evolution equation derived above, 
equation (23). 

On the other hand, the obvious upper limit is furnished by the boiling time calculated in 

To conclude this chapter we would like to point out the region of fulfillment of condition (jH} 
in terms of the heat flux and the size of the channel. Using the definition of the Strouhal number 
(0) and thermodynamical parameters, we obtain that in the temperature region T = 1.4-^2.1K 
condition (jH)) is equivalent to the following one: 

q^L » 0.2 ^ 0.6 -]. 29 

Thus for heat load of moderate intensity q = 1 10 W/cm'^ and for sizes of the setup 
L = 10 -f- lO^cm relation (jHl) and consequently (jSHl) are valid to good accuracy. 

Conclusion 

We have outlined the procedure of separation of the fast from the slow processes in the equa- 
tions of Hydrodynamics of Superfiuid Turbulence. As an illustration of the developed procedure, 
we studied the slow stage of the evolution of the transient heat load of moderate intensity. It 
has been shown how the wave-like behaviour of the heat pulse described by initially hyperbolic 
equations gives way to the parabolic-type non-hnear heat conductivity equation. The latter 
formally coincides with the one derived earlier by Dresner on the phenomenological grounds. 
Unlike the approach advocated by Dressner, our consideration is not only self-consitent but 
also allows to determine a range of validity of the elimination procedure. 
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